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Abstract 

For a system of M coupled Schrodinger equations, the relationship is found 
between the vector-valued norming constants and M + 1 spectra corresponding 
to the same potential matrix but different boundary conditions. Under a special 
choice of particular boundary conditions, this equation for norming vectors has 
a unique solution. The double set of norming vectors and associated spectrum 
of one of the M + 1 boundary value problems uniquely specifies the matrix of 
potentials in the multichannel Schrodinger equation. 
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1 Introduction 

Consider the system of coupled one- dimensional Schrodinger equations 

-^^«(x) + X^K.;3(x)^/3(x) = a = l,...,M. (1) 

In this system, each equation is referred to as a 'channel' and e^s are the energies 
of channel 'thresholds'. Once E > Sa, it is said that a's threshold becomes open. 
The system is a matrix generalization of the ordinary one- dimensional Schrodinger 
equation. The coupled Schrodinger equations originate in the Feshbach's unified theory 
[Ij of nuclear reactions and correspond to so-called approximation of the strong coupling 
(when a finite number of equations in is left). Now, that method, renewed and 
generalized (see, e.g. |2]), finds a lot of applications and, rightfully, is one of the most 
universal tools for microscopic description of systems with many degrees of freedom 
(nuclear structure, reactions, molecules, etc). 

The inverse problem for multichannel Schrodinger equation (^J has also been devel- 
oped [3-5]. As in one-channel case, one can uniquely restore the potential matrix Vajsix) 
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from the spectral measure that, e.g., for the case of bounded interval, is specified by the 
complete set of eigenvalues En and so-called norming vectors (spectral weight vectors) 
7q, (£'„). These vectors characterize the behaviour of the normalized wave functions 
\1/q,(x. En) at one of the boundaries of interval (or at the origin for a half-axis problem, 
etc.), see also below. 

At the same time, in the one-channel case we have more variants of the inverse 
problem. Among them, there is a statement of inverse eigenvalue problem on a bounded 
interval where no norming constants occur. Namely, the potential is uniquely recovered 
from a knowledge of only two different spectra, each for a distinct pair of homogeneous 
boundary conditions (with the same potential) [0] . There were established necessary 
and sufficient conditions of the solvability of the inverse Sturm-Liouville problem from 
two-spectra, see, e.g., the book [Tj. 

Till now, one attempt to generalize this theorem to the multichannel Sturm-Liouville 
operator has been known to the author - the article fS] (the case of a finite-difference 
operator). Though not complete, this work gave an idea of the existence of such a 
generalization in principle. No doubt, the possibility of deriving potential matrix from 
a certain set of spectra would contribute to the multichannel inverse problem theory. 
In present article, results concerning that problem are obtained. It is found that M + 1 
spectra determine the Va/six) and, under special conditions, it is possible to uniquely 
restore multichannel Sturm-Liouville operator. 

The central idea of the paper is to derive the relationship between M + 1 spectra 
and M-component norming vector 'jaiEn) associated with one of the M + 1 boundary 
value problems. Then, having the double set of eigenvalues and norming vectors, one 
can uniquely restore an interaction matrix (by Gel'fand-Levitan procedure). 

Next section is devoted to setting forth these results. We shall find the sought 
expression which, however, does not guarantee the uniqueness in itself. Only under a 
special choice of boundary conditions it can be represented in a form of system of linear 
algebraic equations which give a simple criterion of the uniqueness and solvability. For 
the sake of the reader's convenience, the narrative is organized so that it goes partially 
in parallel with standard derivation of two spectra formulas given in |7 , chapter 3. 

2 Derivation of the formula for norming vector 

We are beginning this section with preliminary notations. Let us rewrite the system 
(P) in a more symbolic form as follows 




(2) 
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where y stands for the whole vector-column solution 



and 



V = Val3 + ea, \ = E. 



The potential matrix is the real symmetric matrix of continuous functions, x G [0, a] 
Next, we add to the equation (j21) the following boundary conditions 



y'{Q) - hy{0) = 
ym - hMO) = 



y'ia) + Hyia)=0 
y[{a) + Hyi{a) = 0, 



M 



(3) 



where we take hi and H to all be the real symmetric matrices. We denote the spectra 
of the M+1 problems © and © by {\n}'^=i and {A^}^]^, respectively. There is no 
theorem of interlacing of the spectra in M-channel case, M > 1. So, we additionally 
require that no spectrum degeneracy should occur. 

Let us denote by A) and Xi{.^i •^) the matrix solutions of the equation Q sat- 
isfying the initial conditions 



0(0, A) = i, A) 1^=0 = Xi(0, A) = i, ^i{x, A)|^=o = h 



(4) 



where the prime stands for the derivative with respect to x. In what follows we shall use 
the prime to denote this derivative apart from the special cases the reader will be let 
know of. Besides, the hat will always stand for the matrix. Do not confuse the following: 
a matrix solution of Q means that each column of the matrix is a vector-solution, only 
satisfying a specific initial (boundary) condition. Eigenvalues of the boundary value 
problems Q and Q coincide with zeros of determinants of the matrices 



1 <l(A) = 0'^(a;,A)U=, + 0(a,A)^^ 

where the bar sign denotes transpose. 

Now we introduce the norming vectors associated with the spectrum {A„} 

/ 71 (A„) \ 



(5) 



oo 
n=l 



7A„ = 



V 7m (An) / 
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so that 



A„)7a„ = y{x, \n) 

with y\x,\n)U=a + Hy{a,Xn) = and Jo" Eali[*a(a;, A„)]2cia; 
spectra {Xi}^^^ 



(6) 

1. Likewise, for the 



/ 7i(AD \ 



7a^, = 



so that 



V 7m (Ai) / 



(7) 



with y'i^x, Xli)\x=a+Hyi{a, X\) = and Y^a=i[^a{x, A^)]^(ix = 1. Let us also introduce 
the function 7a (vs. A) so that 7a = 7a„ when A = A„ and 7a = 7Ajj when A = A^. 
That function makes the sense at the points A„ and AJ^ only. In between, we have the 
freedom to specify it arbitrarily. We can only require this function to be continuous 
different iable and have no singularities. 
We take 



fi{x, A) = 'jxXiix, A) + mi(A)7A0(x, A), 
where mi{X) is scalar and we require that 

/;(x,A)U=, + /,(a,A)^ = =^ 



(8) 



(9) 



mi(A)[7A0'(x, A)|^=a + 7A0(a, X)H] = -[-fxXiix, X)\^=a + 7AXi(a, X)H]. 
Comparing with (jS)) we have 

<li(A)$(A) 



mi(A) 



$(A)$(A) 



where we denote $(A) = $(A)7a and $i(A) = $i(A)7A. 

Next, employing the well known Green formula we have 

/•a ^ no. _ ^ 

(X-Xn) fi{x,X)(t){x,Xnh\^dx = {X - Xn) ^xXi{x, X)(j){x, Xn)l\^dx 

Jo Jo 

<|,.(A)$(A)_ - I ^ 

~(A - An)^7TTir7TT7A0(a:, A)0(x, Xn)i\^dx = f-{x, A)|^=o0(O, A„)7a„ 



$(A)$(A) 



(10) 



-/i(0, A)0(x, A„)|^=o7a„ = Ixik - h)-fx„, (12) 
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where we use the definitions (jH) and (fTTj). The last equahty follows from the fact 
that the matrices in Q are symmetric. 

Let us pass to the limit A — A„. Then the equation (fT^ goes over into 



A — An _ 



$,(A„)$(A„) 



7A„(/i.-/^)7A„ = l, (13) 



where we used the L'Hospital rule. 

We shall prove that this formula can be represented as 

(K - An)-^ U'jh^^xXhr - hhx. = 1, (14) 

where the prime denotes that we omitted, in the product, the term with the number n. 

Since $(A) and ^i{X) are the entire holomorphic functions they are determined (to 
within constant multipliers) by their zeros and, hence, can be represented as follows 

oo \ oo \ 

$(A)=cn(i-i-); $.(A)=an(i-T7)- (15) 

Substituting (fTH|) into (fT!^ we have 



nr=i(i-|f)w 



lxSk-h)lx^ = l. (16) 



Now we have to ascertain the expression for the CC/CiC. We shall need some knowl- 
edge about an asymptotic behaviour of the solutions of (j2I). First of all, these equa- 
tions become uncoupled in the limit A — >■ oo. So, as in one-channel case, we have 
limA^oo$(A){$»(A)}-i = 1, and the same for the transpose of these matrices. Taking 
this into account we obtain 



r^r^ oo \i oo \ \ 

^n^iimn^ = i- (17) 

CiC ^^-^ A^ ^^~;.=i A;, -A 

We have the following asymptotic formulas for A and A': A^ = /i^ + 0(1) and the 
same for A*. Then A^ — A^ = 0(1) and the series X^JT KA^ ~ A^)/(A^ — A)| converges 
uniformly as A ^ oo. Hence, we can pass to the limit in each term of the infinite 
product 

oo \ \ oo \ \i 

•i-n^-iLmn(i+f^)-i. (18) 

We see from and Jllj) that 
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r^r^ oo \i 



(19) 



At last, we can obtain the final expression for 7a„. Substituting (fT^ into (fTBj) we 
have the formula ()14p- the system of M equations {i = 1,...,M) for determining M 
components of 7a„. 

In the one-channel case the formula (jl4j) goes over into the known expression for 
two spectra: 



CO \1 \1 



-(/i2 - /ii)7a„ 



(20) 



where the matrix values become scalars, and we denote, by indeces 1 and 2, two spectra 
determining scalar norming factor 7a„. 

The system ()14j) is not linear one: each row in it contains the quadratic form 7A„(^i — 
/i)7a„. Hence, these equations cannot be solved uniquely in general (including solvability 
itself). In other words, we have to impose some constraint on choosing the matrices hi, 
i.e. the difference hi — h. Among other possibilities, we give several realizations which 
will allow a unique solvability of the system ()14j). 

i) The symmetric matrix hi — h = has the form of a Jacobi matrix: 

\ 



ii) 



11 


S12 





. 





.(i) 

12 





Ai) 
S23 


. 







S23 





Ai) 
S34 





Ai) 

SA/-1A/ 



(21) 



V ... ai-iM 

where the main diagonal contains only one non-zero element, ^[i. Then 



M 



lxSh^ - hh,^ = e}Y7l(An)' + 2 Y: ekUklk-l{Xnhk{Xn). 

fc^l 



(22) 



Introducing the variables ui = 7i(A„)^ and Uk = 7fc-i(Ari)7fc(An), k = 2, ...M we can 
rewrite the last expression as follows 



M 



l\,Xhi - /i)7A„ = ^n^i + 2 ^ iiUk^k- 



(23) 



Then fll4|) becomes the system of linear algebraic equations for the variables uj. If 
t^i = 7i(An)^ > 0, then 7i(A„) = ^uj^^, 72(A„) = ^uj2/ujy^ and so forth. The sign 
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1 /2 

in front of oj^ in the expression for 71 (A„) determines the common sign for 7a„ and, 
hence, is inessential: The whole vector-valued wave function is determined to within 
sign (±). With the non-zero element 7^ 0, / 7^ 1 positioned in arbitrary place of 
the main diagonal, the scheme is analogous. 



ii) The matrix hi — 'h = Cy^ is represented as follows 



/ 



SIX 



V 






■ii) 









Ai) 
SMI 







\ 



Ai) 
SIM 



J 



(24) 



i.e., the matrix contains one non-zero row and one non-zero column which cross each 
other in a place of the entry (li^ . For the quadratic form we have (using the symmetry 
of hi — h) 



M 



(25) 



Introducing new variables 9k = 'yi{\n)lk{^n), k ^ I and 9i = 7i(A„)^ we can now look 
upon (fT^ as a linearized system again: 



M 



(26) 



After deriving 9i, one can obtain 7i(A„) trivially. Of course, the solvability in this case 
depends on whether the corresponding determinant for the system (j26|) is non-zero and 
9i>Q. 

In all the cases, the knowledge of the complete set {A„,7a„}.^i allows a unique 
restoration of the potential matrix by the standard Gel'fand-Levitan theory (its multi- 
channel generalization) . 



3 Conclusions 

In this paper, the relationship is established between components of the norming vector 
7a„ associated with a certain boundary value problem (with the spectrum {A„}J^]^) and 
the spectra (including {A„}J^^) of M -|- 1 multichannel Sturm-Liouville operators with 
the same potential matrix Vapix) but different boundary conditions. As a matter of fact, 
the central result is the formula (fT^ . Though giving no unique solutions in general, it 
can get linear if we require the matrices hi to be of special type. Hence, the uniqueness 
of the multichannel inverse eigenvalue problem from M + 1 spectra is however possible 
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for a particular class of boundary conditions. The problem of specifying the necessary 
and sufficient conditions needs a special examination. It is clear that scrutinizing the 
asymptotic behaviour of the spectra with different boundary conditions will be required. 
It is closely associated with specifying the class of diffcrentiable functions the Vapix) 
pertain to. So, the results given present only an intermediate stage in investigations on 
the subject. 

References 

[1] Feshbach H 1962 Ann.Phys., iVFlQ 287-313. 

[2] Zakhariev B N and Suzko A A 1990 Direct and Inverse Problems (Heidelberg: 
Springer) 

[3] Cox J R 1966 Ann.Phys., Ary39 216-236; 1967 J. Math. Phys. 8 2327-31. 
[4] Cox J R and Garcia H R 1975 J. Math. Phys. 16 1402-09. 

[5] Chadan K and Sabatier P 1989 Inverse Problems in Quantum Scattering Theory 
2nd edn (Heidelberg: Springer) 

[6] Borg G 1946 Acta Math. 78 1-96. 

[7] Levitan B M 1987 Inverse Sturm- Liouville Problems (Zeist, The Netherlands: 
VSP). 

[8] Zakhariev B N 1990 JINR Rapid Communications N6[45]-90 41-48. 



8 



